Dynamical recurrence and the quantum control of coupled oscillators 
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Controllability - the possibility of performing any target dynamics by applying a set of available operations - 
is a fundamental requirement for the practical use of any physical system. For finite-dimensional systems, such 
as spin systems, precise criteria to establish controllability, such as the so-called rank criterion, are well known. 
However most physical systems require a description in terms of an infinite-dimensional Hilbert space whose 
controllability properties are poorly understood. Here, we investigate infinite-dimensional bosonic quantum 
systems - encompassing quantum light, ensembles of bosonic atoms, motional degrees of freedom of ions, 
and nano-mechanical oscillators - governed by quadratic Hamiltonians (such that their evolution is analogous 
to coupled harmonic oscillators). After having highlighted the intimate connection between controllability and 
recurrence in the Hilbert space, we prove that, for coupled oscillators, a simple extra condition has to be fulfilled 
to extend the rank criterion to infinite-dimensional quadratic systems. Further, we present a useful application 
of our finding, by proving indirect controllability of a chain of harmonic oscillators. 



One of the most fundamental questions in science is what 
kind of dynamics a given system can host. Control theory 
addresses this question in the light of how the dynamics of 
the system can change as a response to our attempts of steer- 
ing it. Control theory can be applied at different levels: when 
dealing with computing devices, for example, one could either 
classify their dynamics by the primary logical operations they 
can perform or, at a higher and arguably more useful level, by 
determining what kind of programs they can run. In quantum 
computing the first level is typically determined by the experi- 
ments, and provides one with a description of the Hamiltonian 
of the system under consideration. On the other hand the sec- 
ond level corresponds to the set of quantum algorithms that the 
quantum computer is capable of running. It is at this second 
level that the capability for a device to perform the algorithms 
theorists dream of is established or disproved. 

To connect the experimental and theoretical levels, one 
faces the problem of translating the Hamiltonian description 
to a description in terms of algorithms it can perform. For 
finite-dimensional quantum systems, for instance qubits and 
qudits, this translation has been accomplished in the 1970s 
by the development of an elegant mathematical framework 
dubbed as algebraic control IHQl. In infinite dimension, how- 
ever, such a translation has been so far elusive. Roughly 
speaking, the problem encountered is described as follows 
(see Fig. Q}: in finite-dimensional systems, the state space 
is 'limited', and if the quantum state evolves in one specific 
direction, it will eventually return to where it started. Clas- 
sically this would be Poincare's celebrated recurrence theo- 
rem, whose quantum mechanical counterpart is given in J^. 
Hence, in finite dimension, there is in some sense no need 
to distinguish between opposite directions or, equivalently, 
to specify the direction of time. For control theory this al- 
lows one to develop a picture in which time is fully elimi- 
nated. Single-directed movement in infinite-dimensional sys- 
tems however can carry an 'intrinsic clock', e.g. distance trav- 
elled from some initial state, and time cannot be eliminated. 
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FIG. 1. In finite dimensions, quantum systems recur (a), while, in 
infinite dimensions, single-directed movement can carry an 'intrinsic 
clock' and systems do not recur anymore (b). However, under certain 
conditions, we prove that harmonic oscillator systems (c) can recur 
as if they were finite-dimensional. 

Now in order to assess quantum control of an infinite- 
dimensional system, our question is: do all infinite- 
dimensional systems have an intrinsic clock? A counterex- 
ample is given by harmonic oscillators: just by looking at the 
system state of an oscillator, one cannot tell how long it has 
been running from any specific initial state. Hence, there is 
hope to perform the same time elimination and reach an easy 
description of operations that can be performed in quantum 
harmonic oscillators. Note that these systems are not mere 
theoretical curiosities. The control of the infinite-dimensional 
degrees of freedom of light, trapped particles, nano- and opto- 
mechanical oscillators, superconductors, Bose-Einstein con- 
densates, and of collective spins of atomic vapors or solid- 
state devices are all of major technological interest, and the 
primary way to address most of such degrees of freedom is the 
manipulation of quadratic Hamiltonians, which correspond to 
descriptions in terms of quantum harmonic oscillators 0-01 ■ 
We will see a specific example concerning the control of ar- 
rays of trapped ions at the end of the Letter. 

In this letter, we prove that a restrictive condition has to be 
fulfilled in order to assess the controllability of quantum har- 
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monic oscillator networks. We shall observe that such systems 
share substantial similarities with finite-dimensional dynam- 
ics, which explains the success of previous numerical results 
IToll . We also demonstrate the potential impact of our findings 
by showing how indirect control methods j 1 lj], developed pre- 
viously for spin systems only, can also be applied to oscilla- 
tors, possibly leading to resource efficient cooling and control 
protocols. We start by presenting the basic notions of the al- 
gebraic control, revisiting the proof of the Lie algebra rank 
criterion and finding why it fails to be sufficient for the 
controllability of generic infinite-dimensional systems. Then 
we will focus on quadratic bosonic Hamiltonians, which give 
rise to the so-called Gaussian operations, and we will deter- 
mine a condition such that the rank criterion will still be suffi- 
cient for controllability in the restricted Gaussian sense. In the 
end, we will present an example relevant to arrays of trapped 
ions and to chain of nano-mechanical oscillators, where this 
condition is fulfilled and where local controllability can be 
proven by applying our general analysis. 

Algebraic control - Let us start by reviewing a finite- 
dimensional control setup in quantum physics. Suppose an 
experimentalist succeeds in setting up a system described by 
the Hamiltonian 



H{t) = H + Y / f k (t)H k 



(1) 



k=l 



where the H k are a set of controlling Hamiltonians that can be 
switched on and off. The Schrodinger equation for the time 
evolution operator U then reads 



dU 

~dt 



-iH(t) U U(0) 



(2) 



The main goal of a control theorist is to determine which 
quantum algorithms, i.e. which unitary operators U, the ex- 
perimentalist can, in principle, achieve by setting the right 
switching times for the f k . 

In this Letter we are interested in controlling systems de- 
scribed as coupled oscillators, thus we need to introduce some 
additional notation and terminology. We shall consider an n- 
mode bosonic system, described by n pair of quadrature op- 
erators qj and pj satisfying the canonical commutation rela- 
tion [qk,pi] = iSk.i- By introducing the vector of operators 
R T = (qi,pi, . . . , q n ,p n ), the commutation relation can be 
written as [Rk, Ri] = i£lkl where 57 is the (2n) x (2n) sym- 
plectic form whose matrix elements are Cljk = <5j+i,fc[l — 
(-l) J ']/2 - 6 jt k+i [1 + (-l) J ']/2 in terms of Kronecker deltas 

In particular we will consider systems described and con- 
trolled by Hamiltonians that are bilinear in the quadrature op- 
erators, i.e. that can be written as H = (1/2) i AkiRkRi 
where A are real and symmetric 2n x 2n matrices. The 
corresponding evolution operators in the infinite-dimensional 
Hilbert space, defined as U = e~ lHt , are the so-called Gaus- 
sian unitary operations since they preserve the Gaussian char- 
acter of quantum states Il2l 1 1311 . If we consider the Heisen- 



berg evolution for the quadrature operators vector R, we ob- 
tain the equation U^RU = SR where S = e~ Ant is a 2n x 2n 
matrix belonging to the real symplectic group Sp(2n, It) sat- 
isfying the equation SflS T = O [14j]. Then, restricting to 
quadratic Hamiltonians, there is a one-to-one correspondence 
between the evolution operator in the infinite-dimensional 
Hilbert space U and the finite-dimensional matrix S, in par- 
ticular they correspond to different representations of the real 
symplectic group 01511 . As a consequence Eq. (J2J can be recast 
in terms of the symplectic representation, providing one with 
the following evolution equation 



S(0) = 1, 



(3) 



where, as pointed out above, S is no longer unitary, and 
G(t) = + J2k=i fk(t)AW)Sl is no longer anti- 

hermitian. In this framework, the main goal of the control 
theorist can be summed up by the following question: which 
symplectic time-evolutions S are achievable by controlling 
G(t) via the functions 

Because the solutions of Eq. (01 are elements of matrix 
groups, one can apply the beautiful framework of Lie groups 
to tackle such questions. Let us suppose that, by setting the 
control functions equal to constant values in G(t), we can 
identify a set £ = {G\, • • • , G m } of linearly independent 
generators of a Lie algebra C. This assumption is known 
as the Lie algebra rank criterion. Its relevance is due to the 
fact that, if the corresponding Lie group Q = e c is a sub- 
set of a compact group, then all its elements can be imple- 
mented with arbitrary precision and the system is said to be 
controllable. The Lie algebra rank criterion is an easy and ex- 
tremely powerful criterion for controllability and works well 
for finite-dimensional unitary gates since they are subgroups 
of the compact group SU(n). On the other hand the solutions 
of (0 no longer enjoy this property, so we cannot apply the 
rank criterion directly. 

Why and where is the compactness used to prove that the 
rank criterion is sufficient for controllability? Given a generic 
Lie group Q, a set £ = {Gi, ■ ■ ■ , G m } of linearly independent 
generators of the corresponding Lie algebra, and an element 
K G Q, one can write 



K = e Gltl e G2t2 ■ • -e G "*" 
with Gj £ £ and tj £ M. 



(4) 



In principle, in the expression above, there will be some ex- 
ponentials involving negative times tj, while, if we want K 
to be reachable by control, we need all the times to be pos- 
itive so that every exponential in the product corresponds to 
an evolution described by Eq. ® and obtained by setting the 
control functions equal to certain constant values. As shown 
in JUSQj me compactness of the Lie group is a sufficient 
condition to switch the sign from negative to positive time: let 
G be one of the generators of C and consider a negative time 
t < 0. If Q = e c is compact, then a sequence of positive times 
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tk > always exists such that 



lim 

k— >oo 



„Gtk 



(5) 



In other words, for a given e > and time t < 0, we can al- 
ways find a positive time t > 0, such that, for a given matrix 
norm, ||e GT - e Gt \\ < e. 

One should also notice that this condition is equivalent to say- 
ing that at a certain time, the evolution operator recurs to the 
identity, that is, recurrence is a necessary and sufficient condi- 
tion to revert the sign from negative to positive times. When 
one deals with non-compact groups, the possibility to switch 
from negative to positive times is in general lost. A simple and 
visually clear example in this sense is given by the squeezing 
operation in phase-space |17|] : if we continuously apply the 
squeezing operation, the state gets more and more squeezed 
and recurrence is never achieved. 

However we will show in the following that for coupled har- 
monic oscillators, even considering non-compact Lie groups, 
recurrence takes place with arbitrary precision if an additional 
condition on the system's Hamiltonian is met. Since, as we 
pointed out before, in the proof of the rank criterion, com- 
pactness is used only to revert the sign of time in evolution 
operators, if one is able to achieve this goal by imposing other 
different physical and mathematical constraints, then the rank 
criterion remains necessary and sufficient to prove that the 
reachable set is dense in the Lie group being considered. 

Controllability of quadratic Hamiltonians - Previously we 
showed that, given a system described by Eq. d20l i. a linear 
control problem for the unitary operator U is defined by the 
Schrodinger equation as in Eq. 0. If we consider Hamiltoni- 
ans bilinear in quadrature operators such that 



H k = \Y^ A ^ RsRi k = {°>- 



'}■ 



(6) 



s,t 



the corresponding unitary operators U are infinite- 
dimensional matrices; however we can consider the 
equivalent control problem, with the form of Eq. 0, 
for the finite-dimensional evolution matrix S as 



5(0) = 1 



(7) 



where A(f,t) = A<°> + ££=i f k (t)A {k) - 
The two evolution equations are equivalent and we will focus 
for the moment on the finite-dimensional representation in 
Eq. 10. The symplectic group Sp(2n,R) is a non-compact 
group and thus one cannot apply the rank criterion to assess 
the controllability of the system. Our main result, contained 
in the following theorem, shows that if we can identify a set 
of linearly independent generators of the symplectic algebra 
C = sp(2n, R), such that the corresponding A^ are positive 
definite, one can achieve recurrence with arbitrary precision, 
and thus the rank criterion will still be a sufficient condition 
for controllability. 



Theorem: If A is a positive (negative) definite matrix then, 

Ve > and VT > 0, 3t>T such that \\e~ ASlT - 1|| < e 

where we considered the Euclidean matrix norm ||M|| = 
v/TrfMUT]. 

Proof: If A is a positive definite matrix, because of 
Williamson theorem 112 ill , we can write A = VDV T where 
D = diag-jVi, i>i, . , . , fjv, vn}, Vj £ R + and V belongs to 
Sp(2n, R). This implies that V T VL = VLV^ 1 and thus we ob- 
tain An = vDnv- 1 . 

The matrix DVl is a normal matrix diagonalized by a uni- 
tary matrix U, such that Dfl = UD'W with D' = 
diag{+ifi, —iv\, . . . , +ivN, —Wn}- Then we have 



An = vud'u^v- 1 = WD'W~ 



(8) 



that is Afl has pure imaginary eigenvalues and is diagonalized 
by the matrix W = VU. As a consequence, we can write the 
matrix S(t) = e ~ Ant as 



S(t) = WE{t)W- x 



(9) 



with E(t) = diag{e- Wl ,e Wl , 
to 



, e 11/11 , e Wn }, which leads 



||s(*)-i|| = || WW) -i)^- 1 !! (io) 

< \\W\\ WW-'W \\E(t)- 1||. (11) 

The matrix W is independent on time t, and thus 
||W|| HW^ 1 !! = K is constant. Let us consider the remaining 
term 



\E(t) in = ^n\E(t) ip] 

= (2^|e-^*-l| 



1/2 



(12) 
(13) 



fc=i 



It is easy to check that \\E(t) — 1|| is a sum of trigonomet- 
ric exponential functions and thus is a quasi-periodic func- 
tion. Because of this property, a time r > T exists, such that 
||-E(t) — 1|| is close to zero with arbitrary precision 12211 . In 
particular, for a given e > 0, we can choose a time r such that 
||£(t)-1|| < e/A' and then obtain the thesis || S(r)-l|| < e. 

m 

It is worth also to notice that our result can be extended to 
semi-definite matrices A's such that QA is diagonalisable, but 
it cannot be extended to generic semi-definite A's. For in- 
stance, the Hamiltonian H = p 2 for a single degree of free- 
dom would not recur. 

This theorem assures that if, by properly choosing the con- 
trol functions ff. (t) in the Hamiltonian given in Eq. d20l i, we 
can identify a set of linearly independent generators of the 
symplectic group, such that the corresponding matrices A^ 
are positive (negative) definite, then we can always revert the 
sign of negative times in the expression corresponding to Eq. 
(0. As a consequence the Lie algebra rank criterion remains 
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a necessary and sufficient condition to asses the controllabil- 
ity of the symplectic group, even if the group is not compact, 
and thus can be used to assess which Gaussian operations can 
be realized given a certain control problem as indicated in Eq. 
©. 

On a more fundamental level, our argument highlights a hith- 
erto unnoticed connection between the normal mode decom- 
position of positive definite quadratic Hamiltonians - formally 
an implication of the Williamson theorem - and the controlla- 
bility of sets of coupled oscillators. This connection is bridged 
by the notion of dynamical recurrence which, regardless of the 
infinite-dimensionality of the Hilbert space, is always guaran- 
teed for positive definite quadratic Hamiltonians. The general- 
ity of our result makes it at first easy to overlook its potential 
impact: the applicability of quantum control to continuous- 
variable systems paves the way to vastly improving fideli- 
ties of current experiments as well as using control more ef- 
ficiently. This fact can be clearly seen in the simple but quite 
surprising example we provide below, where our result is used 
to simplify the controllability properties of a large harmonic 
oscillator network. 

Local controllability of a quadratic harmonic oscillator 
chain - One of the most important requirements in quantum 
information and quantum computation is to dynamically ad- 
dress and control individual interacting systems. From an ex- 
perimental point of view, it is also desirable to have complete 
control on a large network by acting only on a small part of 
it. Indirect control has been already proved for qubit systems 
1 1 111 , while, regarding networks of harmonic oscillators, it has 
been recently shown, for instance, that by probing only one 
site of the network, one can reconstruct the full quantum state 
of the system 11811 . Here we use our theorem to prove the con- 
trollability of a chain of harmonic oscillators, where only one 
or few sites of the chain are accessible. This kind of example 
is relevant if we think, for example, of an array of interacting 
trapped ions, where in principle one can implement Gaussian 
operations by addressing single ions and manipulating their 
trapping frequencies |Q uj^ . 

Let us start by defining the bosonic mode operators dj = 
(qj + ipj)/\/2 and a] = (qj — ipj)/\/2, satisfying the com- 
mutation relation [a$, aj] = 8ij. Let us consider an n-mode 
bosonic chain, described by a Hamiltonian as in Eq. d20l i. In 
particular the always-on Hamiltonian reads 



H 



,t. 



n-l 

92 ^(ajaj+i + h.c). 

3 = 1 



n-l 



a 3 a ]+l 



h.C. 



(14) 



where, for the sake of simplicity, we consider all the oscil- 
lators having with the same frequency to. If we consider 
.9i = 92, this corresponds to the qjqj+i coupling, which is the 
most common between the harmonic oscillators' interactions. 
From now on we will consider the renormalized coupling con- 
stants cjj = gj/cu and assume that they are both positive; a 



sufficient condition for the positivity of Hq (for every number 
of bosons in the chain n) is gi + g 2 < 1/2. We consider as the 
controlling Hamiltonians, a local phase-rotation and a local 
squeezing term on the first mode of the chain only, i.e.: 



Hi 



(jJi[a\aA and H 2 — x( a f + a \ 2 ) 



(15) 



We prove that, by denoting with C the symplectic algebra, 

C = (iH ,iH 1 ,iH 2 ) [ .,. ] , (16) 

that is, by computing all possible commutators of these oper- 
ators, of any order, and their linear combinations, we can ob- 
tain all the elements of C (details of the proof can be found in 
the Supplemental Material). Since the Lie algebra is a vector 
space, any set of linearly independent linear combinations of 
the above operators satisfies the rank criterion. We have then 
to show that, by properly setting the control functions fk{t), 
we can identify one of these sets, such that the corresponding 
matrices are positive definite. There are in principle infi- 
nite choices, in particular it is easy to check that the following 
set fulfils all the conditions above 



Hq — Hq 

H 1= H + aH x 

H 2 = H + pHi + 5H 2 



(17) 

auji > (18) 
< S X < Pwi ■ (19) 



In practice, this example is directly relevant to arrays of 
trapped ions and chains of nano-mechanical oscillators. For 
instance, in the case of the transverse ionic modes, local con- 
trols analogous to ( |22] | could be obtained by manipulating the 
local trapping frequencies, as detailed in 11911 and realised in 
Jfl (analogous forms of control have been envisaged for opto- 
mechanical setups as well I20I1 ). Note also that even restricted 
local control might suffice for certain manipulations, depend- 
ing on the desired tasks. For instance, as a side product of the 
proof reported in the Supplemental Material, one can show 
that, if the local control is restricted to phase rotations gener- 
ated by Hi, the whole symplectic algebra is not achievable, 
but all the passive operations, comprising beam-splitters and 
local phase rotations can be realized. This would allow one to 
implement, for example, cooling protocols based on swapping 
excitations between sites of the array. 

Conclusions - While the general control theory of infinite- 
dimensional systems remains hard, we have found a surpris- 
ingly simple solution for the case of quadratic interactions of 
coupled harmonic oscillators. To demonstrate the applicabil- 
ity of our result, we have also discussed its application to indi- 
rect control (and, potentially, cooling) of chains of oscillators. 
It is worth mentioning that proving controllability and find- 
ing an actual control pulse are completely distinct tasks. For 
instance, using the quantum recurrence theorem on the the- 
ory level is useful, but for a control pulse one could not rely 
on it, as the recurrence and hence the resulting pulses would 
take far too much time. This is well understood and typically 
overcome by using numerical routines to optimize the pulses. 



In our case, a similar point arises regarding the requirement of 
positive (negative) definiteness of the Hamiltonian, which was 
used as a sufficient element of the proof. In an actual pulse se- 
quence, it could be beneficial to use non-positive or negative 
definite Hamiltonians in order to achieve faster control. 
Last, but not least, let us remark that the dynamics of clas- 
sical systems governed by quadratic Hamiltonians is also de- 
scribed by the symplectic group of canonical transformations: 
our finding hence applies, as it stands, to the controllability of 
classical, as well as quantum harmonic oscillators. 
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SUPPLEMENTAL MATERIAL 
Proof of rank criterion for harmonic oscillator chain 

Here we consider the control problem defined by the Hamil- 
tonian 



H(t)=H + ^2f k (t)H k , 



(20) 



k=l 

where the always-on Hamiltonian reads 

n / * \ n—1 

H =ojY^ \ a ) a 3 + 2 ) + 51 + hx -)- 

3=1 ^ ' 3=1 

n-1 

+ 92 y^{aj(ij + i + h.c), 

3 = 1 



(21) 



the local controlling Hamiltonians are 

i?i = wi (a\ ai \ and H 2 = x(af + a\ 2 ) , (22) 



and dj and a] are n pairs of bosonic operators satisfying the 
commutation relation [aj,afc] = 5j. k - 

In the following we prove that the Lie algebra rank criterion 
for this system is fulfilled, that is, by denoting with C the sym- 
plectic algebra, then 



C = (iHo, iHi,iH; 



2/h-]- 



(23) 



We will work with the infinite-dimensional representation of 
the real symplectic group, where the basis of the correspond- 
ing Lie algebra reads 



C ={i{aja\ + ajafe), a 3 a\ ~ aja^, , 
i(ajc4 + aja k ), a\a\ - dja k } 
with j, k = 1, . . 



(24) 



Proving the Lie algebra rank criterion corresponds to show- 
ing that by computing all possible commutators of the control 
operators H k , of any order and their linear combinations, we 
can generate all the elements of the basis of the algebra listed 
above. For the sake of simplicity we will prove it by consid- 
ering the case gi = g 2 in the Hamiltonian Hq d2Tb ; this corre- 
sponds to considering a qjqj+i coupling between the sites of 



the chain, while the most general case (which comprises also 
the rotating-wave approximation case where <?2 = ) can 
be proved following the same line of reasoning. We start by 
showing that we can generate all the elements corresponding 
to the two first sites of the chain: 

= \[iH 2 ,iH 1 ] = af -a\ (25) 

0\ = [iHo, iHi] = a\a) 2 + a\a 2 — a\a\ — a\a 2 (26) 
2 = [iHi, 0{\ = i{a\a\ + a\a 2 + a\a\ + 0102) (27) 

R ( a 2) = \ [([O 1 ,iH ]+2O 2 ),iH 1 ] = a\a 2 - a x a\ (28) 

(29) 
(30) 



T^=O x -R^=a\a\-a 1 a 2 



SW = [T^\R^]+SU = 



,t2 



Tb = i°i - R { a 2) ^ H i\ = M4 + aia 2 ) (31) 
3 = [iHx,Oi -2i?i 12 >] 

= i(a\a\ — a\a 2 — a\a\ + ciia 2 ) (32) 

Rb 2) =\(O a - Os) = i(4a 2 + oioj) (33) 

pW = \(KH l + [R£ 2 \RW])= l a\a 2 (34) 

(35) 



2 \-— ' 1 b 
S^ = [P^\S^]=i{af + al) 



The operators P< 2 ), s£\ S?\ S?\ T< 12) , if 2) , R ( ^\ 

(12) 

R b , together with the local control operators iH k , complete 
all the operators corresponding to the first two sites. Because 
of the symmetry of Ho, we can proceed in the same way for 
all the operators belonging to the neighbouring sites. Then, to 
obtain the two-mode long-distance operators, one can easily 
show that 



[i{a' k ai + a k aj), ajaj. - ajafc] = i(a]a/ + a a\) (36) 

and analogue commutators for the remaining terms. The proof 
is hence complete. ■ 
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